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Total  progeny  in  a critical  age-dependent 
branching  process  with  immigration 


by  Howard  Weiner 


1.  Introduction. 

At  time  t •=  0,  a renewal  process  starts  with  I.I.D.  interarrival 

times  having  non- lattice  distribution  function  G (t),  G (0+)  ■ 0.  These 

o o 

epochs  are  the  arrival  times  of  new-born  immigrating  cells,  where  k cells 
arrive  with  probability  Pj^^,  and  let,  for  0 < « < 1» 


(l.l) 


h (a)  - S p.  8 , 
® IfO  “ 


and  denote 


(1.2) 


0 < 0 - h'(l)  - 2 kp.  < 
k-1 


Each  iosnigratlng  cell,  independent  of  any  other  cells,  initiates  an 
age-dependent  branching  process  [3]  with  cell  lifetime  distribution  G(t), 

G(0f)  ■ 0,  and  non-lattice.  The  offspring  generating  function  is 

m 

h(s)  * ^ p.s  . Assume  that  each  initiated  branching  process  is  critical, 
lo»0  “ 

that  is. 


(1.3) 


h’(l)  • 2 kp.  ■ 1. 


Dmot*  by 


(1.4)  Z(t)  “ the  total  number  of  cells  born  by  t arising  from  all 

cells  Immigrating  by  t and  their  respective  Initiated  critical 
age-dependent  branching  processes. 

The  purpose  of  this  paper  Is  to  obtain  an  explicit  limit  law  for 
the  Laplace  transform  of  Z(t)  by  using  the  corresponding  result  for  the 
Galton-Watson  or  discrete  time  process  obtained  by  Fakes  and  a series 
of  approximations. 

2.  Integral  Equations  and  Approximations. 

Assume  that  at  t * 0.  there  are  no  cells  present  due  to  limnlgratlon, 
and  hence  no  new  cells  to  Initiate  a branching  process. 

Ut 


(2.1)  N(t)  " total  nuntber  of  cells  born  by  t In  a critical  age- 
dependent  branching  process  as  given  by  (1.3) > (1.4). 


Define 


(2.2) 


1(0, t)  ■ E exp-(0N(t)). 


Then 

(2.3)  f(0,t)  - a“*(l  - C(t)  + j\(#(0,t-u))dC(uy]  . 
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(2.4) 


F(®,t)  - E exp- (©2(0), 


where  2(t)  is  es  In  (1.4). 

Then  f or  t > 0 , arguing  as  in  14], 

(2.5)  F(G,t)  » 1 - G^(t)  + J^h^($(e,t-u))F(©,t-u)dG^(u) 

and  F(O,0)  "I. 

Let 

(2.6)  ^ + J^h(«^(©,t-u))dG(u)],  n 

and  «g(0,t)  ■ 1. 

Lee 

*(©,irt-l)  - e'’®h(S(0,n))  for  n > 1 


(2.7) 


L«t  a > 0 b«  • conatanc.  Daflna  tba  approxlmanta , for  n>  I, 


(2.8)  - 1 - 

and  ^ 

Similarly.  l«c 

(2.9)  p(©,(n+l))  - h^(»(e,[na]))P(e,n)  for  n > 1, 

and  P(8»0)  = 1, 

which  implies 

n 

P(0,nfl)  - n h (i(©,(ma))). 

BSil  ® 

Ws  note  the  relationships  ((6],  lensm  3 and  also  p.  285) 

(2.10)  *(«.t)  i *(8)  as  t - • (see  also  [51,  p.  228) 

(2.11)  #(«.n)  4 *(®>  as  n -*  •, 

where  #(©)  is  the  transform  of  a bona  fide  random  variable. 

(2.12)  n(l  - *(®/n^))  < K < •.  all  n, 

•iMre  K la  a eemstaat. 


U 


3.  Limit  Theorema . 


Theorem  1 . Assume  that 


(3.1) 


0 < h(0)  < 1,  h^(0)  < 1, 


and  that 


h'(l)  - 1,  0 < h"(l)  < «,  h^(l)  » 

Denote  y -'|h"(l),  and  a = p/y. 

Then  for  all  0 > 0,  a > 0, 

(3*2)  11m  P(-~,n)  = (sech  a ^y9  . 

Proof . We  will  Indicate  the  adaptation  of  the  method  of  proof  of 
Pakes  [6]  applied  to  this  more  general  model,  as  Fakes'  proof  applies  only 
to  the  case  where  the  lasalgratlon  mean  Interarrlval  time  - mean  lifetime  - 1 
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' ’ r —i. 


The  nodal  here  corresponds  to  Che  case  where  the  lomlgration 
lacerarrlval  Cine  is  a > 0,  and  the  mean  lifetime  is  1. 


Case  I.  * > !• 

Using  Fakes'  notation  ([6}  eq.  20, p. 285)  for  out  case. 


(3.3) 


P(0,(iri-1))  - n h^(«(e,tma])). 

ja»l 


Since  <[6]  eq.  21,  p.  285) 


*(0»k)  #(0)  as  k “*  •, 


we  may  write,  via  ([6]  eq.  30,  31,  p.  285),  where  0 = ©/n 

n 


(3.4) 


log  P(©„,(tH-l))  - -PE  (l-«(e„,tmal)+R^"\e) 

® BS«1  ” 


where  R'  '^(0)  ”♦  0 as  n • if  E (l>f(0  ,[ma])  is  bounded  in  n. 

« n 


Writing 


(3.5) 


-P  E (l-4<©  ,[ma])  » -p  E (l-6(©  )) 
B^l  " n-1  “ 


+ P E (♦(«„JnaJ)-*(0  )), 

V*1  ” ^ 


and  by  Fakes  ([6]  eq.  31,  p.  287), 


(3.6) 


-P  E (l-«(«  ))  --WO/Y)^'^. 


Ualng  tht  proof  of  L«na  4,  pp.  266>267  and  p.  288  of  [6],  mm  obtalaa 
Ch«c,  aa  n • 


(3.7)  8 Z («(e„,(ma])-«(e^)) 


~ii  S -i, 

Y n , , , -6rma]/n  V J ^ . 


Y n ^ Jq  ^^-6ax 


where»  following  [6],  we  have  denoted 
(3.8)  6^2  . 


(3.9) 


•1  \ r 1+e"**  > 


eq.  (3.4)  • (3.9)  yield  the  result  for  the  case  a > 1. 
jCesjs^^II.  0 ^ a ^ 1. 

Observe  first  that  if  l/a  is  a positive  integer,  that 


(3.10) 


-[na]  -,1/a 

P(G,nfl)  - n h (t(9,k)) 

^k-1  ° 


If  a ■ f/r,  where  0 < i < r are  relatively  prime  positive  integers, 
then  it  can  be  seen  that 


(3.11) 


.'ln*l  N 

P(9,tH-l)  -V  nh^(*(9,k))j 


In/r] 

n h (»(9,ki)). 
iMl  ° 


Now  let  n « in  (3.10),  (3.11),  replacing  9 by  9/n^.  Apply  Pakes' 
liait  result  ([6],  p.  285)  to  the  rl^t  side  of  (3.10),  and  to  tha  first 
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product  on  the  right  of  (3.11).  Apply  the  result  of  Case  I above  to 
the  second  product  on  the  right  of  (3.11).  This  proves  Case  II  for 
0 < a < 1,  a rational. 

Suppose  0 < a < 1 is  arbitrary.  Choose  a sequence  of  rationale 
(a^}  such  that  a^  **  a. 

By  the  (uniform)  continuity  and  monotonicity  in  a of  the  product 
for  P(0,n+1)  given  on  the  right  side  of  (2.9),  it  follows  that  one  can 
choose  the  convergent  sequence  of  rationale  [a^}  such  that  for  all 
r > R^»  n > N^,  either  (by  a slight  expansion  of  notation) 

(3.12)  P(0/n^,n,a^)  < P(©/n^,n,a)  < P(e/n^,n,a^j) 


or 

P(e/n^,n,a^_j_^)  < P(©/n^,n,a)  < P(©/n^,n,a^) . 

Now  let  n -•  *,  then  r • in  (3.12),  iising  the  continuity  of  the 
function  on  the  right  side  of  (2.2).  This  suffices  for  the  result. 

Lenma.  Let  G^(t)  and  G(t)  two  distribution  functions  such  that 

G^,(Of)  - G(Of)  - 0 

and 


0 < ■ I tdC(t)  < •, 

1 •'o 


Assume 


0 < 1 t^dG(t)  < « 

*'0 

0 < 1 t^dG  (t)  < ». 

••o  ° 

Let  and  G^(t)  denote  the  r-th  convolutions  of  G^  and  G,  respectively. 

Let  G * H(t)  denote  the  convolution  of  G and  U.  Then,  uniformly  In  k, 

0 < k < n, 


(3.13) 

11m  G, . , 

t 

* G , (t)  = 0 

o,n-k'  ^ 

when 

for  all  € > 0, 

(3.14) 

lim  fl  - 
t -»  0# 

when 

for  all  e > 0. 

Proof.  Let  {X^}  be  I.I.D.  each  with  distribution  G. 

Let  {Y^}  be  I.I.D.  each  with  distribution  G^,  and  Independent  of 
Define 


(3.15) 


S 


n,k 


<l^] 

L 


(X^-EXp  + 


n-k 


Then  to  prove  (3.13),  (3.14)  together,  lu  suffices  to  show  that 
(see  [2]) 

(3.16)  llm  P(ls„  . I > 3cnu)  - 0, 

uniformly  in  k,  0 < k < n. 
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By  Chebyshev's  inequality. 


(3.17)  < 


[kalVar  4-  (n-k)Var  X + n Var  Y 


(3em^n)‘ 


(3€ni2n) 


which  goes  to  zero  Independent  of  k,  as  n -*  <». 

Note:  By  using  the  techniques  of  [1],  it  may  be  possible  to  relax  the 

second  moment  conditions  and  still  obtain  that  Ptjs  . 1 > e]  goes  to  zero 

n I K 

uniformly  in  k,  as  n -•  »,  but  this  has  not  been  done. 


Theorem  2.  Assume  the  conditions  of  Theorem  1 and  the  lemma.  Then 


amj^/m^ 

lim  F(-^,t)  * fsech 
t - - t^  “l  -* 


Proof.  The  proof  is  by  use  of  approximants,  writing 
(3.18)  F(e,t)  - P(9,(n+D)  s F(9,t)  - + ^'(nfl)^®’'^^  " P(9,(n+1)). 

<3-19)  F^^^j(9,t)  - F(9,t)  - Jjh^(«[^3(e.t-u))F^(9.t-u) 

- h (i(0,t-u))F  (9,t-u)  + h ($(0,t-u))F^(9,t-u) 

o ft  o n 

- h^(S(Q,t-u))F(0,t-u)]dG  (u). 

o o 


n- 
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llie  proof  will  be  divided  Into  several  claims,  using  an  abbreviated 
notation.  Only  the  proofs  of  Claims  IV,  VIH  will  be  given,  since  the 
others  either  follow  similarly,  or  are  similar  to  those  given  in  [2],  [6] 


(3.20) 

.a. 

Claim  I: 

n 

(3.21) 

Claim  TI: 

^(n+l) 

(3.22) 

Claim  III; 

(3.23) 

Claim  IV; 

"n+1  ■ 

C . , 1 ★ G ^ 1 I ( 

[ka)  o,nfl-k 


+ G (t). 
o,n+l  ' 


Proof.  By  (3.19), 


Vl  - " ■ K <'’o<*lna])  - »o<‘» 


Applying  the  mean  value  theorem  to  the  first  term  on  the  right  of  (3.24), 
bounding  by  1 in  the  integrand  yields 

(3.25)  + . G^(t> • (!-*(«)) . 

Iteration  of  (3.25)  and  using  Fj^-F-1-F<lin  the  last  step  yields 
the  result. 


(3.26)  Claim  V;  *j^(©,t)  - «(0,n)  > 0,  all  n,  0,  t, 

(3.27)  Claim  VI;  F„(©,t)  - P(©,n)  > 0,  all  n,  ©,  t. 


M.2.)  CUlsHI:  *„,-•(««<  (1-0  .n„>o. 


(3.29)  Claim  VIII; 


P^l(0,t)  - P(g.tH-l)  < c(l-i(©)).\\i_G,^  ,)  * G ft) 


where  c = maxfl.m^)  and  G_^  = 0 for  0. 


Proof . By  the  expression 


+ . h^($([naJ))]dG 


+ (1-G  (t))(l-P(n)), 


«nd  the  mean  value  theorem,  (3.30)  yields 


(3 


Vl  " < J(V*’("»<lGo  + J 


^*[na]  - *([n*]))dG^  + C 


By  (3.28)  and  (3.31), 


»•“>  Vi  - + ..a  - . o„(,,  , 

U.r.tlo„  «ul  t|„  . 0 » . ,(2)  < 

yield  the  result.  ~ 


l-C„(t». 

.(«)• 

1 - e.(t) 
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We  may  now  complete  the  proof  of  Theorem  2. 
Claims  IV,  VIII  yield  the  Inequalities 


< P(«.c)  - P(e,n+1)  < c(i-»(e))  z a-o,..,)  » G„  „ . 
- - [kaj  o,n-k 


+ c(l-G  ) * G 

' o'  o,n-l 


r t(l+c)  -| 

Now,  set  ""m J tight  (left)  inequalities,  respectively. 


Substitute  6/n  for  6 in  (3.33).  Then,  since  from  (2.12), 


1 - *(0/n^)  < ^ 


for  some  constant  K,  let  t ->  «<>.  The  lenana,  via  uniformity  of  approach  of 
the  summands  to  zero,  yields  that  the  right  and  left  sides  of  (3.31)  go  to 
zero.  Now  Theorem  1 completes  the  proof. 
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